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1 .  Introduction 

This  report  describes  the  design,  construction  and  evaluation  of  a  working 
torsion  disc  viscometer. 

Initial  tests  at  gas  pressures  between  0.1  mbar  and  1000  mbar  show  that  the 
theory  of  the  instrument  is  adequate  and  that  the  viscometer  is  capable  of 
detecting  small  changes  in  gas  viscosity. 

A  thorough  investigation  of  the  effect  of  radioactive  irradiation  on  gas 
viscosity  will  be  made  shortly,  following  the  work  of  Kestin  and  Shah. 

An  analysis  of  their  approach  is  given  here,  viewed  in  the  light  of  our 
preliminary  work,  and  suggests  the  most  likely  ways  of  achieving  a  change  in 
gas  viscosity. 

2.  Review  of  Work  to  Date 

Kestin  and  Shah7(see  also  Clark,  Kestin  and  Shah1)  made  an  experimental 

investigation  of  the  changes,  caused  by  ionisation,  in  the  apparent  viscosities 

of  various  gases.  They  used  a  torsion  disc  viscometer  to  measure  the  viscosity 

and  radioactive  sources  to  produce  ionisation  of  the  gas.  The  changes  in 

apparent  viscosity  were  generally  small,  the  most  significant  being  decreases 

of  0.25  %  and  7.6  %  for  air  at  1  atmosphere  and  1  mm  Hg  pressure  respectively, 

and  an  increase  of  5.25%  for  argon  at  1  atmosphere. 

To  investigate  the  possibility  of  changing  viscosity  under  flight  conditions, 

we  constructed  an  aerodynamic  skin  friction  drag  balance,  which  could  measure 

the  drag  on  a  flat  plate  for  air  flows  of  up  to  200  m/s  at  atmospheric  pressure. 

High  accuracy  (y  0.1%)  was  achieved  at  nigh  speeds,  but  no  changes  in  drag 

on  irradiation  were  observed.  At  low  speeds  the  accuracy  is  limited  due  to  the 

small  drag  forces  (down  to  0.1  gf)  and  these  errors  were  greater  than  or 

2  4  15 

comparable  with  the  magnitude  of  the  change  expected  ’  ’ 

Accordingly  we  decided  to  make  our  own  measurements  on  almost  static 
gases  using  a  viscometer  giving  higher  accuracies  than  the  drag  balance. 


2. 


The  type  of  viscometer  chosen  was  the  torsion  disc  viscometer  with  a  free  disc 
suspended  from  a  fine  torsion  wire  (figure  1).  This  is  simpler  than  rotating 
cylinder  or  fixed  disc  types,  and  is  easier  to  use  when  making  an  extended 
series  of  measurements  under  different  conditions. 

It  was  also  the  same  choice  as  made  by  Kestin  and  Shah,  although  our  approach 
makes  the  use  of  the  viscometer  rather  more  straightforward. 


3.  Theory 

Consider  a  circular  disc  of  finite  but  small  thickness  (i.e.  h  <<  R) , 
immersed  in  a  fluid  of  density  p  and  viscosity  y  and  suspended  from  a  thin 
torsion  wire  giving  a  restoring  force  X0. 

If  the  disc  is  set  oscillating  the  amplitude  will  be  damped  by  the  viscous 
action  of  the  fluid  (the  wire  will  also  contribute  to  the  damping,  this  is 
considered  later). 

The  equation  of  motion  of  the  disc  can  be  written 


16  +  D  (0,0)  +  X0  =  0 

where  the  damping  is  some  function  of  0,9  and  the  properties  of  the  fluid. 
If  the  v.  ping  is  small  (D(0)<<  X)  and  roughly  a0  (fluid  velocity  small)  the 
motion  will  approximate  to  a  damped  simple  harmonic  motion. 


Writing 


10  +  D0  +  X0  =  0 


the  damping  can  be  characterised  by  the  logarithmic  decrement  of  the  system 


A  =  In 


(-  )- 

\ 0n+ 1  / 


where 


An  experimental  determination  of  A  will  therefore  give  us  information  on 
the  fluid  provided  we  can  find  an  expression  for  D. 


3. 


(An  expression  for  T  can  easily  be  found  in  terms  of  the  properties  of 
the  suspension  wire  - 

g 

According  to  Landau  and  Lifshitz  the  restoring  couple  for  a  rod  twisted  about 
the  axis  of  symmetry  parallel  to  its  length  is 


ce 

£ 

The  shear  modulus  S 

S  = 


where  C  =  J| 

is  related  to 
E 

2(l+o) 


Sirr 


(3) 


the  Young's  Modulus  and  Poisson  ratio  t/ 


(A) 


and  this  leads  to 


T  =  16 
o 


/  7T  (  1+0  ) 

V—Ed*5" 


+o)Il 


for  the  natural  period. 

The  period  with  damping  ,  T,  is  given  by 


T  =  T  /I  +  (£  )z  -  T  (1  +  {  &)2) 
o  2ir  o  zir 


(5) 


(6) 


for  simple  harmonic  motion) . 

To  find  an  expression  for  D  we  need  to  solve  the  equations  of  motion  of 
the  fluid  around  the  disc.  The  Navier  Stokes  equations  take  the  following  form 
in  cylindrical  polar  coordinates  (neglecting  pressure  gradients  and  assuming 
that  the  motion  is  cylindrically  symmetric,  valid  assumptions  for  small  damping 
after  the  decay  of  any  transients). 


(7a) 


(7b) 


(7c) 
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(7d) 


The  presence  of  the  non  linear  terms  make  a  solution  difficult.  However, 
to  get  an  idea  of  the  type  of  motion  let  us  suppose  that  Ur  and  Uz  are  small 
and  that 

D  +  0  as  z,r  +  » 


U^(z=±— ,r<R)  =  rw0ocosu)t  =  l^coswt 


i.e.  we  neglect  the  damping  intially;  6q  is  the  amplitude  of  the  motion. 


The  equations  then  reduce  to 


au* 

p— — <P  =  y 

■at 


32U 


y4>  +  - 
3rz  r  3r 
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Near  the  surface  of  the  disc 


U,  «  ruicoswt 


and  so 


!?“♦  +  i  !“♦  .  "4  .  o 


3r 
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suggesting  that  in  the  region  of  the  disc  a  reasonable  approximation  to  the 
equation  of  motion  is 


314  32U, 

c  H*  ‘ " 


or 


is*  .  vi^ 

at  3^ 


(9) 


This  linear  equation  can  be  readily  solved  by  the  separation  of  variables 
.  13 

(equivalent  to  Stokes  2nd  problem  ):-  restricting  our  attention  to  the  top 
of  the  disc,  taking  z=o  on  the  surface  for  simplicity,  we  put  U  a  f(t)g(z)  to  get 


1  df 


1  u  C  . 

,  =  v — *r  =  SI  where  SI  is  some  complex  constant. 

f  dt  g  dZ 


Sit 


M 

=  e»  v 


neglecting  the  arbitrary  multiplicative  constants 


Then  1=6,  g 
for  the  moment. 

To  fit  the  boundary  conditions  at  the  disc  we  put  SI  =  iw 


Jl  - 


then 


(1+i) 


5. 


Rejecting  the  +ve  root  to  match  the  boundaiy  condition  at03  we  get 

.  „  (l+i)2 

,,  lwt  2v 
V  “  e  e 


or,  taking  the  real  part 


-z/X 


V  =  rwGcos  (cot  -  z/X)e 
where  X  =  /— -  =  I-—  is  the  penetration  depth. 

V  U)  v  irp 


(10) 


The  solution  takes  the  form  of  an  oscillation,  damped  in  the  z  direction, 
so  that  the  region  of  interest  is  within  a  few  X  of  the  disc. 

The  surface  shear  stress  at  the  disc,  is  given  by 


T  = 

v 


z=o 


=  (sin  oi t  -  cos  mt) 

A 


00 


(»e-e) 


The  term  in  0  can  be  thought  of  as  effectively  increasing  the  inertia  of 
the  disc,  due  to  the  fluid  carred  by  the  disc.  For  small  damping,  this  can 
be  neglected  and  our  expression  for  the  total  damping  is 


D  =  2 


'o 


p  0  3. 

Y  2it  r  dr 
A 


(2  sides) 


Including  a  correction  for  the  edge,  assuming  a  similar  equation  of  motion, 
we  get 

y_uoip 


D  =  itR't  / ( 1  +  ^L) 


and  the  logarithmic  decrement  is  thus 


A  =  ttR4/ttpPT  (1  +  2h) 
21  R 


(12) 


The  damping  due  to  the  wire,  A  ,  can  simpiy  be  added  to  this  expression  to 
get  the  total  damping,  provided  Aq  is  small. 

This  expression  is  in  agreement  with  measured  drags  for  our  viscometer 
to  within  a  factor  of  1.5  or  so,  underestimating  the  actual  damping. 

The  major  part  of  the  discrepancy  can  be  attributed  to  the  effect  of  the 


edge  of  the  disc. 


A  better  expression  has  been  obtained  by  Maries  &  van  Paeniel 


Writing 


6 


10 


U  =  r  \  (r,z,t)  where  \  is  the  angular  amplitude  of  the  fluid,  they  solve  the 

T 

sinplified  Navier  Stokes  equation. 


with  the  boundary  conditions. 

'jz(t)  =  0  (t)  =  0(O)e^t  on  the  surface  of  the  disc  and  'fCz.r  ■*■«>)  =0. 

The  time  dependence  is  assumed  to  be  as  for  damped  SUM,  i.e.  ft  =  iu)  -  A/T. 

The  solution  for  Che  damping  can  be  pieced  together  by  considering  different 
regions  of  the  fluid  in  turn. 

Near  the  edge  (region  !  of  figure  1)  it  can  be  assumed  there  is  no  z  dependen 
and  using  the  time  dependence  this  gives 


+  2  dil 


r  dr 


+  p2  ^j=  0 


where  p2  - - . 

r  v 

writing  \  =  f (r)  we  see  that 
r 


d2f  +  J_  df 
dr2  r  dr 


+  (p2 


(14) 


or,  with  p  =  pr 


d  f 
dp* 


4—  +  ( t  -  J2) 

dp  pz 


f  =  0 


which  is  Bessel's  equation  of  order  1. 

The  solution  for  4"i  can  therefore  be  expressed  as  a  linear  combination  of 
first  order  Hankel  functions  of  the  first  and  second  kind:- 


CiHi(l)(pr) 

r 


+  C2H(2^ (Pr) 
r 


L 


t 


% 


7. 


To  meet  the  bondary  conditions  at  00  and  on  the  disc  we  need 

\  =  0(t)  R  H{l)(pr)  (15) 

rHi ' ^  (P r) 

Off  the  corner  of  the  disc  we  can  separate  variables  and  use  \,=  \\  for  z  =  h/2  to  get 

\  =  Q(t)elp(2~h/2)RHr1)  (pr)  06) 

rHi'^pR) 

Above  the  disc,  assuming  no  r  dependence,  but  allowing  for  the  influence  of  the 
fluid  in  region  2,  the  equation  of  motion  is 

I  ii3  =  lib*  i  1^2 

v  3t  Sz"  R  3r 

9 

With  the  identity 

sailisi;:) ,  -  iilOo  ,  Ho) 

dr  r  1  o 

this  gives  a  differential  equation  in  z  which  is  readily  solved  to  give 

\  =  0 (t)  {1  -  C 2 i  (z-h/2)/R]  [  |  -  H^1)(pR)/H1(1)(pR)]  }  eip(Z_h/2) 

3  pr  o 

Knowing  \  ,  the  viscous  torque  on  the  disc  can  be  calculated,  3nd  after  some 
algebra ,  expanding  the  Hankel  functions  in  their  asymptotic  series,  Mari ens  and 
Van  Paemel  get  the  following  expression  for  the  logarithmic  decrement:- 


r=R 


A  =  IL'gmPt  [1  +  f  +  A  (2  +  |i)  +  |  (A)2 


0  +  ~)  *  0((A)4)]  (19) 


For  small  X/R  it  can  be  seen  that  this  reduces  to  our  earlier  expression. 

The  viscosity  can  be  obtained  from  A  by  interation  or  by  using  a  look-up 
chart . 

.10  .  1 

Experimental  data  obtained  by  Marions  and  van  Paemel  and  by  Kestm  and  Wang 

show  that  this  equation  generally  gives  values  of  y  correct  to  a  few  percent, 

entirely  satisfactory  for  measurements  of  small  viscosity  changes. 


8. 


A  more  accurate  theoretical  expression  would  undoubtedly  be  difficult  Lo 

obtain,  as  centrifugal  force  must  cause  a  radial  motion  of  the  fluid  near  the 

disc  and  consequent  axial  flow,  i.e.  and  will  not  be  negligible. 

An  approximate  solution  for  the  problem  of  a  constantly  rotating  disc  with 

1  ?  11 

U^,  U_^  4  0  has  been  obtained  (see  Sparrow  and  Cregg  “  and  Schliehting  )  but 
not  for  the  case  of  an  oscillatory  motion.  However,  adapting  the  results  to  an 
oscillating  disc  suggests  a  similar  magnitude  for  A. 

'V 

For  A  >  R  our  expression  for  A  must  break  down  due  to  the  asymptotic  expansion 
of  the  Hankel  functions.  This  will  happen  for  lowT  gas  pressures.  An  approximation 
to  the  solution  can  be  obtained  by  solving  equation  9  in  a  similar  fashion  to 
before,  but  with  u  =  o  on  the  walls  of  the  viscometer. 

This  gives 

A  =  [  C  (b/\)+  C  <V\)+  ~  C  (y.\)  ]  (20a) 


where  C.(b)  =  cosh  (b/A)  sinh(b/A)  +  cos  (b/A)  sin  (b/A) 
sinh2(b/A)  +  sin2  (b/A) 


For  large  b/A  this  reduces  to  equation  12  and  for  smell  bA  reduces  to 


A  = 


ttR4Tp  +  1  +  4h  L 
A I  bj  b2  R  "b  3 


(20b) 


This  gives  the  correct  dependance  on  y,  although  the  discrepancy  between 
the  magnitude  of  the  measured  drag  and  A  calculated  from  20b  is  rather  worse  than 
at  high  pressure.  This  is  probably  because  the  assumptions  involved  in  equation  9 
are  unrealistic  when  A  'v  R;  in  particular  the  influence  of  the  side  walls  of 
the  viscometer  is  likely  to  be  of  great  importance.  An  analysis  of  this  problem 
will  be  made  soon.  However,  the  use  of  an  equation  of  the  form 
A  =  Ky 


where  K  is  an  empirically  determined  constant  characteristic  of  the  disc  should 
be  adequate  for  measuring  changes  in  apparent  viscosity  at  low  pressure. 


9. 


The  approach  used  by  Kostin  and  Shah  in  their  work  also  involved  an  attempt 
to  solve  t He  simplified  Navier-Slokes  equation.  Rather  than  assume  the  form  of 
the  time  dependence,  the  Laplace  transform  was  used  to  replace  the  time-dependence 
of  the  Navier  Stokes  equation  and  of  the  equation  of  motion  of  the  disc.  The 
viscous  torque  is  constructed  assuming  that  the  torque  is  proportional  to  the 
torque  calculated  assuming  the  solutions  for  an  infinite  disc  and  an  infinite 
cylinder.  The  "constant"  of  proportionality  is  an  empirically  determined  function 
of  the  bondary  layer  thickness. 

The  resultant  inverse  Laplace  transform  cannot  be  solved  exactly,  but  after 


some  analysis  the  expression 


c((£>2>  ♦  i  (i  ♦  f)  J-  (i  -  |a  -  Ia2)’-^5  c(|) 

l  K/2  K/Q  2BIK 


- 2  <!  -  v  -  ° 

where  0  =  T/Tq  and  <5  -  /v/w  was  obtained  by  Kestin  and  Wang^  and  was  used  together 
with  an  empirically  determined  correction  function  by  Kestin  and  Shah^  to 
calculate  apparent  viscosities. 

This  proceedure  warrants  some  criticism;  firstly  the  viscosity  only  enters 
the  problem  through  the  correction  function,  which  has  values  up  to  3  for 
measurements  at  atomospheric  pressure.  This  makes  it  hard  to  justify  the  use 
of  such  a  complicated  expression,  with  all  its  computational  difficulties.  Secondly, 
the  equation  does  not  approach  the  correct  limit  when  6/R  is  small. 

We  therefore  feel  that  using  the  Mariens  -van  Paemel  equation  (19)  and 
equation  (20b)  for  A>R  is  a  better  approach.  A  possible  further  theoretical 
solution  might  be  achieved  by  assuming  the  time  dependence  of  the  solution  and 
Fourier  or  Laplace  transforming  wrt  z.  However,  the  inverse  transform  is 
again  the  difficult  step  (see  ref.  3  for  a  general  discussion  of  the  problem). 


A .  Design 

Having  chosen  the  type  of  viscometer,  the  choice  of  disc  and  torsion  wire 
material  is  critical.  Initial  studies  led  to  a  disc  of  several  cm.  diameter,  with 


a  thickness  of  about  Irani.  The  wires  used  were  restricted  to  those  readily 
available,  being  either  steel  (piano  wire)  or  CuBe,  of  'v  0.2  ram  diameter. 

Neither  material  showed  any  significant  deterioration  with  time,  and 
the  restoring  forces  were  accurately  linear.  The  wire  length  used  was  generally 
between  10  and  20  cm.  Some  experiments  with  quartz  fibres  were  made,  but  due 
to  their  brittleness  these  were  not  used.  There  were  also  doubts  about  their 
long  term  strength.  This  could  be  affected  by  water  vapeur  for  example.  The 
chosen  wire  was  fastened  to  the  disc,  and  also  suspended, by  means  of  pinchucks. 

To  measure  the  angular  position  of  the  disc  a  small  galvanometer  mirror  was 
fixed  to  the  base  of  the  lower  pinchuck,  just  above  the  disc.  The  beam  from 
a  Spectra  Physics  1  mW  He/Ne  laser  was  reflected  from  the  mirror  to  a  screen 
'v  2  m  distant.  The  mirror  was  5  mm  in  diameter  with  a  focal  length  of  1  m. 

This  gives  a  better  spot  on  the  screen  than  with  a  plane  mirror;  about  1  to  2  mm 
across  typically.  No  sizeable  amplitude  errors  are  introduced  by  the  curvature  of 
the  mirror  provided  the  beam  is  central  and  the  mirror  is  not  too  far  from  the 
axis  of  the  rotation. 


Starting  the  torsional  motion  in  a  consistent  manner  and  without  causing 
excessive  non-axial  vibration  was  a  problem.  This  has  been  overcome  by  using 
an  automatic  starting  system.  This  consists  of  a  d.c.  electric  motor  acting  on 
the  upper  pinchuck.  The  motion  is  opposed  by  a  spring,  with  a  rigid  stop  defining 
an  accurate  zero  position.  A  circuit  has  been  built  to  give  a  preset  number  of 
pulses  of  the  motor  at  the  resonant  frequency  of  the  viscometer  (figure  2).  This 
gives  a  consistent  starting  amplitude  and  causes  minimal  other  disturbance  of  the 
disc. 

The  disc  is  housed  in  a  ^  6  mm  thick  glass  cylinder  of  18.3  cm  internal 
diamter  and  32  cm  length.  The  ends  have  ground  glass  flanges  which  seal  by  "0" 
rings  to  1  cm  thick  dural  plates.  The  upper  plate  has  gas  light  lead  tli roughs 
for  a  thermocouple  and  the  starter  motor.  A  hole  in  the  plate  connects  to  a 
Schaeritz  PTD  3I0W  pressure  transducer,  an  Edwards  Model  17  Pirani  gauge  and 


11. 


the  gas  /  vacuum  lines.  The  lower  plate  is  held  by  spring  clips  and  is  easily 
removed  for  the  introduction  of  radioactive  sources  (figure  3).  These  are 
placed  on  the  base  and  irradiate  the  underside  of  the  disc.  There  is  no  need 
to  alter  the  disc;  this  eliminates  any  change  in  calibration.  The  laser,  screen 
and  viscometer  can  be  aligned  accurately  using  the  various  reflected  beams  from 
the  glass  cylinder  and  galvanometer  mirror. 

The  viscometer  is  contained  in  an  insulati  ng  box  having  thick  polystyrene 
sides  and  a  polymethylmethacrylate  (PMMA)  front.  The  temperature  can  be  controlled 
thermostatically  (using  light  bulbs  for  heating)  to  an  accuracy  of  0.05°C. 

To  achieve  maximum  accuracy  and  ease  of  operation,  we  decided  to  measure 
the  logarithmic  decrement  (log.  dec.)  automatically.  Instead  of  measuring  a 
series  of  amplitudes  by  eye,  a  circuit  was  built  to  count  the  number  of 
oscillations,  n,  needed  for  the  motion  to  decay  to  various  known  amplitudes. 

The  log.  dec.  can  then  be  obtained  from  the  slope  of  a  graph  of  JhiO  versus  n. 
Provided  n  (which  is,  of  course,  only  accurate  to  the  nearest  integer)  is  large 
enough,  the  accuracy  is  good.  By  making  four  counts  from  an  initial  known 
amplitude,  a  five  point  graph  is  obtained.  The  log.  dec.  used  is  such  that  there 
are  between  70  and  130  counts  per  section  in  normal  use,  giving  an  accuracy 
of  'v  0. 1  %.  As  this  is  comparable  with  the  natural  variation,  this  procedure 
does  not  cause  any  loss  in  accuracy.  The  experimental  configuration  is  shown 
in  figure  A.  The  phototransistors  (PT's)  are  used  to  detect  if  the  amplitude 
on  the  screen  is  greater  than  the  PT  position.  Although  the  transistors  have 
lenses  of  3  ram  diamter  they  trigger  if  the  edge  of  the  laser  spot  reaches  the 
edge  of  the  lens.  The  accuracy  in  amplitude  is  thus  limited  to  the  accuracy 
in  measurement  of  the  PT  position,  typically  ±  0.23  mm  to  ±  0.5  mm.  The  positions 
of  the  phototransistors  were  chosen  to  give  roughly  equal  intervals  in  JcnO. 

6  is  calculated  from 

tan  (0+60)  =  ^ 
d 

,  o 

The  screeen  is  perpendicular  to  the  laser  beam  to  <  1  ,  a  negligible  error. 

The  zero  error  is  calculated  from  the  time  spent  by  the  laser  beam  to  the  left 
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The  centre  phototransistor  can  also  be  used  for  measuring  the  period  of  the 
oscillation. 

Figure  5  and  6  are  schematic  diagrams  of  the  circuit  used.  Figure  7  and  8 
show  the  details  of  the  major  sections  of  the  circuits  (i.e.  details  of  displays 
etc.  are  omitted).  An  IC  list  and  pinouts  are  given  in  the  appendix. 

The  operation  of  the  circuit  is  as  follows:  the  PT's  are  connected  to  a  trigger 
(555)  which  produces  a  voltage  pulse  of  a  few  tenths  of  a  second  duration.  This 
gives  reliable  triggering  and  also  drives  an  LED  for  a  visual  check  on  the 
triggering.  The  central  PT  drives  a  3  digit  BCD  counter  via  a  divide  by  two 
(JK  flip  flop).  PTO  resets  the  counter,  and  PT  1  to  PT  4  update  the  corresponding 
displays  to  the  current  state  of  the  counter.  Thns  each  display  is  updated 
until  the  laser  beam  just  fails  to  reach  the  PT,  and  the  count  is  then  frozen. 

A  "master  latch  enable"  (MLE)  output  is  derived  from  PT4  to  freeze  the  display 
a  couple  of  minutes  after  the  last  pulse  from  PT  4  (i.e.  at  the  end  of  a  run). 

This  prevents  any  spurious  signals  causing  a  loss  of  data  if  the  system  is  left 
for  some  time.  The  start  button  puts  MLE  high  and  also  triggers  the  starter 
circuit. 

The  timing  circuitry  operates  in  one  of  two  modes,  zero  position  or  period. 

In  the  zero  position  mode  the  time  spent  to  the  RHS  and  LHS  of  PTC  are  measured 
by  enabling  two  6  decade  counters  for  an  appropriate  interval.  The  counters  are 
clocked  at  10kHz(nom.)  by  a  555,  giving  timing  to  0.1  ms. 


In  period  mode  succesive  periods  are  measured  on  alternate  displays, 
each  counter  being  reset  after  the  display  has  been  held  for  half  a  cycle. 
An  output  is  available  for  calibration  of  the  clock  frequency,  if  required. 
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In  zero  position  mode  the  counters  are  reset  by  one  of  PT 0  to  PT  4,  to  give 
a  measurement  at  a  particular  amplitude. 

The  counters  drive  multiplexed  displays  via  buffer  transistors;  the  strobe 
output  is  also  used  to  strobe  the  common  cathode  of  the  displayed  counts  to  give 
similar  brightness. 

Various  sources  of  error  must  be  considered  in  measuring  the  log.  dec.  in 
this  way.  Firstly,  consider  errors  in  the  calculated  zero  position.  We  have 


assumed  constant,  whereas  it  will  change  by  'v  over  one  oscillation.  Since 


60  =  0cos 


T  +T 
L  R 


•  •  I  .  w  . 

this  gives  'v  5%  error  in  60,  or  6x.  As  Sx  <  5  mm  typically, 


this  error  is  negligible. 


An  error  in  the  timing  of  and  TRwill  arise  from  the  finite  size  of  the 

phototransistor;  the  timing  sequence  is  indicated  in  figure  9.  If  x  -  a  sin  tut 

r~5 - 2  ^  dv  -  <ux  tux  .  . 

v  =  atu  cos  out  =  uj/a^-  x*1  =>  -r~  =  ~2_'"Z  ~ - z  at  x  as  x  <<  a.  The  timing 

QX  3  “03  3  Z  Z 

error  arises  from  the  incorrect  travel  time  across  half  of  the  pt,  due  to  the 

1  dx  dx 

difference  in  velocity  over  half  the  pt:  i.e.  dt  'v  -r  (— —  -  -  ) 

4  V  V+dv' 

1  dxdv 

<v  - - 2_ 

4  \ 

'v  dx;  xvdx  (23) 

4  a  vz 

'v  1  x7dx2 
4  a^tu 


and  with 


x  ^5  mm 


dx  'v  3  mm 


w  'v  1 


a  'v  200  mm 


we  get  dt  ^  1  ps 

which  is  quite  negligible. 

The  effect  of  the  finite  laser  spot  size  should  be  of  a  similar  magnitude. 


(The  value  for  6x  generally  varies  by  a  few  tenths  of  a  mm  between  runs, 

which  can  probably  be  attributed  to  small  timing  errors  due  to  slight  non-axial 

vibration  of  the  disc).  Secondly,  a  small  correction  due  to  the  PMMA  can  be 

made  to  the  amplitudes  x  (figure  10) 

m 

6  =  t  sine  (1-jjj)  (24) 

where  N  is  the  refractive  index  of  PMMA  =  1,495.  6  ranges  up  to  ^  0.7  ran  for 

the  m  =  o  position.  Thirdly,  the  beam  path  on  the  screen  is  roughly  parabolic 
due  to  the  vertical  inclination  of  the  incoming  beam. 

The  beam  height  at  the  screen  is 

H  *  (1  +  Kl)2)H  (25) 

d  o 

where  Hq  is  the  height  of  the  beam  at  the  central  pt.  This  could  give  errors 
due  to  triggering  at  positions  other  than  the  edge  of  the  PT.  However,  if 
is  minimised  the  maximum  deviation,  H-Ho,  is  only  a  few  mm,  and  this  can  be 
compensated  for  by  a  small  vertical  displacement  of  the  PT.  Finally,  the  main 
cause  of  error,  external  vibration,  must  be  considered.  Although  this  cannot 
create  a  torque  about  the  wire  axis ,  if  the  non-axial  vibration  is  sufficient 
to  cause  a  few  mm  displacement  at  the  screen  there  is  a  chance  of  slight 
inaccuracies  in  the  counting.  If  the  displacements  are  greater,  then  there  is 
the  possibility  of  missing  a  count.  This  is  less  important  for  PT  0-PT  4  as 
they  only  need  1  in  2  "attempts"  to  suceed,  and  provided  the  final  contact  is 
made,  i.e.  when  amplitude  =  position  of  PT,  no  error  is  introduced.  For  the 
centre  PT,  the  high  velocity  of  the  spot  seems  to  give  better  reliability.  On  the 

f  the  screen  is  set  up  correctly  it  is  possible  to  keep  vertical  vibration 
amplitudes  £  1  mm,  and  the  circuit  usually  performs  with  100%  reliability. 

Measurements  at  low  pressure  can  be  more  difficult,  as  the  log.  dec.  for  pendulum 
oscillations  is  also  reduced.  However,  it  is  obvious  from  the  results  when 
a  significant  number  of  counts  have  been  missed,  and  the  log.  dec.  can  be  found 


15. 


5.  Experimental  Results 

For  each  experiment  the  four  numbers  and  two  times  recorded  by  the  automatic 
circuitry  are  fed  into  a  Hewlett  Packard  HP41CV  programmable  calculator. 

This  calculates,  and  corrects  for,  the  zero  position  error  and  then  perrorms 
a  least  square  regression  of  the  counts  (including  count  0)  on  the  natural 
logarithm  of  the  angular  amplitudes. 

Both  A  and  the  goodness  of  fit  are  calculated.  Sample  results  are  given 
below  (zero  position  calculated  using  PT  2) 


36  Gauge  Piano  Wire 

counts  96  192  289  387 

T,  and  T„  2.3255  2.3810 

L  K 


38  Gauge  Cu/Be  Wire 
74  149  223  297 

5.5289  5.6082 


The  parameters  calculated  from  these  figures  are:- 

6  =  -6.04  mm  6  =  -3.62  mm 

x  x 

A  =  0.004099  A  =  0.005341 

r2  =  0.999985  r2  =  0.999996 

(see  figures  11  and  12) 

In  a  series  of  several  measurements  A  can  typically  be  determined  to  ~  0.1  %> 
this  is  of  the  desired  order  of  accuracy. 

A  series  of  measurements  on  air  at  various  pressures  and  at  20.05  ±  0.05°C 
has  been  performed. 

The  specifications  of  the  wire  and  disc  used  are  as  follows: 

Disc 

m  =  47.15  g  I  =  4.28  x  10  ^  kgm^ 

R  =  42.49  mm  (not  corrected  for  pinchuck; 

h  =  1.00  mm  not  an  appreciable  error), 

centre  hole  3  mm  radius  (OBA  clearance) 


Wire 

Cu  /Be 

d  =  0.148  ±  0.001  mm  E  =  1.17  x  1011  Pa  (ref.  14) 

1=16  cm  a=0. 35  (ref. 14) 


The  vi scone  ter  was  set  up  with 
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bi 
b2 
b  3 

d 

X 

o 

XI 
x2 
X3 
Xl* 


1820  ±  10  mm 
750.0  mm 

488.5  mm 
323.0  mm 

214.5  mm 
143.0  mm 


to  nearest  0.5  cm 


to  nearest 
corrected 


0 . 5  mm; 
for  perspex 


J 


In  leak  measurements  the  viscometer  was  found  to  leak  from  vacuum  at  a 
rate  of  'u  0.02  mbar/hour.  The  measured  period  was  11.27  s.  Equation  5  leads  to 
a  value  of  11.5  ±  0.3  s.  The  results  of  the  tests  are  presented  in  figure  13 
with  the  wire  damping,  Aq,  taken  as  0.00045  ±  0.00005. 

It  can  be  seen  that  there  is  very  good  agreement  between  experiment  and 
equation  19,  within  the  region  of  validity  of  the  equation.  At  lower  pressures 
equation  20  gives  the  correct  viscosity  dependence.  Below  1  mbar  the  logarithmic 
decrement  appears  to  fall  again  although  the  results  showed  some  irreproducibility. 
This  could  be  due  to  slip,  effusion  of  lighter  gases  into  the  system  or  an 
actual  change  in  viscosity.  This  is  to  be  investigated;  a  slip  or  viscosity 
change  seems  unlikely  as  the  mean  free  path  is  no  more  than  a- few  mm  (see  Kennard'’)  . 

Generally,  the  results  indicate  that  the  theory  of  the  viscometer  is  adequate 
and  that  it  can  be  used  to  measure  viscous  drags  over  a  wide  range  of  gas 
pressure . 


6 .  Comments  on  the  Experimental  Results  of  Kostin  and  Sh a h 

The  results  of  Kcstin  and  Shah  generally  showed  few  significant  changes 
in  viscosity.  The  most  significant  changes  were  observed  at  low  pressure,  with 
a  disc  that  was  calibrated,  removed  and  coated  with 


P0  (an  a  emitter)  and 
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replaced  in  the  viscometer.  Slight  increases  in  damping  were  observed  for  most 
gases,  but  a  reduction  of  7.6  %  in  apparent  viscosity  was  observed  for  air  at 

1  mmHg. 

As  at  most  a  few  measurements  were  made  for  each  gas,  the  magnitude  of 
the  error  was  uncertain,  but  probably  was  of  a  few  %.  This  means  the  increases 
are  not  at  a  high  confidence  level.  There  is  also  the  possibility  of  a  change 
in  calibration  having  occured.  Another  problem  is  that  the  pressure  varied 
by  up  to  several  percent  between  the  radioactive  and  inactive  measurements. 

For  example,  there  was  a  5%  decrease  in  the  density  of  the  air  for  the  1  mm  Hg 
measurements.  Since  the  viscometer  may  have  been  operating  in  a  pressure 
dependent  regime  (cf  figure  13)  this  might  have  been  a  significant  error. 

In  the  other  case  of  significance.  Argon  at  1  atmosphere,  a  +5.2%  change 
was  observed,  with  a  1%  increase  in  pressure. 

7.  Conclusions  and  Future  Work 

A  very  accurate  viscometer  has  been  developed  and  built  to  measure  small 
viscosity  changes  over  a  range  of  pressures.  Our  initial  measurements  together 
with  an  analysis  of  Kestin  and  Shah's  work  suggest  that  air  and  argon  should 
be  studied  over  a  range  of  pressures,  the  low  pressure  measurements  being  conducted 
in  the  pressure  independent  regime.  A  more  accurate  theory  of  the  low  pressure 
behaviour  should  also  be  developed.  To  investigate  possible  chemical  effects 
on  irradiation  measurements  should  be  made  with  inert  gases  and  with  gas  mixtures 
containing  potentially  reactive  gases  such  as  ethylene.  An  investigation  of  the 
effect  of  water  vapour  concentration  would  also  be  of  interest. 
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Appendix 

IC  List  and  Pinouts 


TIL  313 

7  Segment  LED  display 

555 

Timer 

4001 

Quad  2i/p  NOR 

4002 

Dual  4 i / p  NOR 

401 1 

Quad  2i/p  NAND 

4016 

Quad  Switch 

4027 

Dual  JK  Flip  Flop 

4069 

Hex  NOT 

4070 

Quad  2i/p  EX-OR 

4510 

BCD  Counter 

4511 

BCR  7  Segment  Display  Driver 

50395 

Six  Decade  Counter/Display  Driver 

4001  BE 


4002BE 


401 1  BE 


1 


4069U8E 


4070BE 


TIL  313  DISPLAY 


Schematic  Diagram 
Viscometer  . 


Fig.  5.  Schematic  Diagram  of  Data -Logging 
Circuit  Counters. 


Fig.  6.  Schematic  Diagram  of  Dat 


Circuit  •—  Timers  . 
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Fig. 9.  Zero  Pc 
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Fig. 10.  Correr 
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11.  Typical 
Steel  Tor: 
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